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I. INTRODUCTION 

In the second-order phase transitions, the correlation 
length £ diverges, and as a result various thermodynamic 
quantities diverge [l[ . In the dynamic case (dynamic crit- 
ical phenomena) , the relaxation time of the order param- 
eter also diverges, which is known as the critical slow- 
ing down. In the context of the AdS/CFT duality 0- 
H| , the dynamic critical phenomena have been studied in 
Refs. [6r{lll. 

In the dynamic critical phenomena, singular behaviors 
are not limited only to the quantities related to the or- 
der parameter. In general, if the system has conserved 
quantities (such as the energy- momentum tensor T**"), 
the dynamic universality class and the dynamic critical 
exponent change due to the coupling between conserved 
charges and the order parameter (mode- mode coupling). 
As a result, one has singular behaviors in the other trans- 
port coefficients (such as the shear viscosity 77Q) 

To be more specific, let us specify a class of theories 
we consider. Below we consider a class of gauge theo- 
ries with local SU(N c )x global U(l) symmetry. So, the 
gauge theory has a conserved charge associated with the 
global C/(l) symmetry, and we consider the case where 
the charge is an order parameter of a second-order phase 
transition. As we discuss below, this class includes the 
Af = 4 super- Yang-Mills (SYM) theory at a finite chemi- 
cal potential and QCD. According to the classification of 
Hohenberg and Halperin [l2j (see Sec. Ill Cl for the classi- 
fication), the dynamic universality class of these theories 
is model H. Model H has a conserved order parameter 
with additional conserved quantities such as T^ w . Model 
H is a typical universality class for a field theory since 
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1 In this paper, we use fj for the shear viscosity to avoid confusion 
with a static critical exponent n. 



a field theory has the energy-momentum tensor. This 
class shows singular behaviors in the charge conductiv- 
ity A and the shear viscosity fj due to the mode-mode 
coupling (Sec. MI Bp . 

We point out that this mode-mode coupling is sup- 
pressed by 1/N% for large- N c gauge theories. Namely, 
this effect cannot be seen in the large- N c limit, i.e., the 
gauge coupling gym — > 0, N c — > 00 with a fixed large 
A := g YM N c . From the gravity point of view, this mode- 
mode coupling effect cannot be seen in the classical su- 
pergravity approximation, and it should emerge from a 
one-loop computation. We illustrate this point in detail 
using the example of the Af = 4 SYM. 

The Af = 4 SYM can undergo a second-order phase 
transition at a finite chemical potential. The gravity dual 
of the system is known as the five-dimensional R-chargcd 
black hole (which is also known as "spinning" D3-brane 
or STU black hole) [TH-EHl- Its dynamic critical phe- 
nomenon has been studied in Refs. 0, ■ In this system, 
the charge density shows a singular behavior, and the dy- 
namic universality class is model B. For model B, only the 
order parameter shows a singular behavior (Sec. Ill Dl) . 
but the system obviously has T^ v . Thus, one expects 
that the universality class becomes model H, but this 
is not the case in the AdS/CFT duality (see Sec. IIII A|) . 
This is because the leading-order AdS/CFT computation 
cannot see the mode-mode coupling effect. 

In this paper, we mostly focus on the Af = 4 plasma, 
but our argument is generic and is independent of the 
details of gauge theories. In Sec. IIVBI we discuss the 
implications to QCD. Son and Stephanov have argued 
that QCD belongs to model H at the critical end point 
Since the model H effect is suppressed by l/N%, we 
discuss the implications to the heavy-ion collisions. 
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II. DYNAMIC UNIVERSALITY CLASS 
A. Hydrodynamic variables 

In the study of critical phenomena, one begins to iden- 
tify the order parameter and to study the statics. How- 
ever, in order to determine the dynamic universality 
class, the full set of macroscopic variables (hydrodynamic 
variables) may become important as we will see below. 

Typical macroscopic variables are 

1. conserved charges such as charge, energy, and mo- 
mentum densities 

2. Nambu-Goldstone modes (if there is continuous 
symmetry breaking) 

3. order parameter (which arises via critical slowing 
down in a second-order phase transition) 

These modes may not be all independent: some con- 
served charge may become an order parameter. Also, 
not all components of and T^ u correspond to hydro- 
dynamic variables. Only conserved charges are hydrody- 
namic variables: they are guaranteed to survive in the 
hydrodynamic limit uj(q) — > as q — > because of con- 
servation laws. To isolate hydrodynamic variables, it is 
useful to carry out the tensor decomposition of J M and 
T^ v . For example, J M can be decomposed as the longi- 
tudinal diffusive mode and the transverse modes. The 
charge density p appears in the diffusive mode. Simi- 
larly, for T^ v 1 the energy and longitudinal momentum 
density appear in the sound modes and the transverse 
momentum densities appear in the shear modes. 



B. Static universality class 

In static critical phenomena, various static quantities 
associated with the order parameter diverge at the crit- 
ical point. The divergences are parametrized by static 
critical exponents. These exponents define a static uni- 
versality class. 

Traditionally, there are six static critical exponents. 
For ferromagnets, they are defined as follows: the spe- 
cific heat: Oh oc le^l - ", the spontaneous magnetiza- 
tion: m oc \ex\^ (T < T c ), the magnetic susceptibility: 
X oc |£t| -7 , the critical isotherm: m oc \h\ x / s {T = T c ), 
the correlation function (T ^ T c ): G{r) oc e~ r ^, the 
correlation function (T = T c ): G(r) oc r - d »+ 2 ->) j an d the 
correlation length: £ oc |eT| _ ". Here, m is the magnetiza- 
tion, h is the external magnetic field, £t " (T — T c )/T c , 
and d s denotes the number of spatial dimensions. For the 
Ginzburg-Landau theory, 

{a,p,iAv,ri)= (0,^,1,3,^,0) . (2.1) 



These exponents are not all independent, but they satisfy 
static scaling relations: 

a + 2/3 + 7 = 2 ,7 = 0(5-1) ,7 = v{2--q) ,2- a = vd s 

(2.2) 

except the last hyperscaling relation, which often fails 
(Sec. IIV C| . The scaling relations suggest that there is 
some structure behind them, which is known as the scal- 
ing law. 



C. Dynamic universality class 

In dynamic critical phenomena, the relaxation time r 
of the order parameter also diverges as r oc £ z , which is 
known as the critical slowing down. The exponent z is the 
dynamic critical exponent. The effect of slow dynamics is 
not limited to the order parameter. The other transport 
coefficients also have singular behaviors in general due to 
couplings with the order parameter. 

The details of the dynamic exponent depend on dy- 
namic universality classes. The dynamic univers ality 
classes were classified by Hohenberg and Halperin (l2j |: 
they are known as model A, B, C, (D, E), F, G, H, and 
J0 We review model A and B in the next subsection and 
review model H in Sec. IIII Bl The classification is based 
on 

1. whether the order parameter is conserved or not, 

2. whether there are the other conserved charges, 

3. whether there exists couplings between these 
modes. 

The dynamic universality class depends on additional 
properties of the system which do not affect the static 
universality class. In particular, conservation laws play 
an important role to determine the dynamic universality 
class. A conservation law forces the relaxation to pro- 
ceed more slowly. As a consequence, even if two systems 
belong to the same static universality class, they may not 
belong to the same dynamic universality class. 



D. Models with order parameter only 

Let us start with the simplest systems where only the 
order parameter matters near the critical point. These 
are model A and B. 

For model A, the order parameter is not conserved 
(such as the uniaxial antifcrromagnet). In this case, the 



2 Model D and E arc written in parentheses since model D reduces 
to model B and model E is a special case of model F with an 
additional discrete symmetry. 
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equation of motion for the order parameter takes the form 
(uj: frequency) 



-iTo 



(2.3) 



The relaxation rate Yq — > at the critical point because 
the Hamiltonian has a flat potential there. Thus, the 
relaxation time t ~ l/To diverges at the critical point. 
For the Ginzburg-Landau Hamiltonian, To oc (T — T c ) oc 
£~ 2 , so t oc £ 2 , namely z = 20 

For model B, a conserved charge, say p, is the order pa- 
rameter (such as the uniaxial ferromagnct). A conserved 
charge satisfies the diffusion equation (q: wave number) 



uj = —iDq 



(2.4) 



where D is the diffusion constant. The diffusion constant 
D is related to the susceptibility \ :— dp/ dp (p: chemical 
potential for the charge) and the conductivity A as 



D = X X 



(2.5) 



This equation has the following interpretation. On one 
hand, the charge fluctuation p diffuses by the charge dif- 
fusion J = — Z?Vp. The same process can be regarded 
as the charge conduction under an electric field E: 



-DVp 



-D 



dp 
dp 



V/i = XE 



(2.6) 



Namely, the charge fluctuation p{x) generates a chemical 
potential p(x) which gives E = — V/i. From Eq. (|2.6p . 
one gets Eq. (|2.5p . 

For model B, A remains finite at the critical point. 
Because x diverges at the critical point, D — > and r ~ 
-(Imw) -1 ~ (Dq 2 )^ 1 diverges. In order to extract the 
dynamic critical exponent, set g£ ~ 1 in the dispersion 
relation (|2.4|) . and determine z from lu ~ £ _z . Using the 
definition of static exponents 7, v and a scaling relation 
"1, one gets x ~ V jv ~ t 2 ~ T ' '• Thus, 



lu = —i\x <z 



r 2+ Y 



-(4-7,) 



r (4 ~ r,) (for qt ~ 1) , 



so 



'/ 



(2.7) 
(2.8) 

(2.9) 



The logic behind setting q£ ~ 1 (scaling form) is as 
follows. We use hydrodynamic dispersion relations such 
as Eq. (|2.4[) . It applies in the hydrodynamic regime 
(q£ <C 1) and takes into account the only lowest power 
of q. However, in the critical regime (q£ ^> 1), higher 
powers of q in the dispersion relation arc no longer neg- 
ligible. In this regime, write u> ~ q z . (The correlation 



3 More generally, when the system has the anomalous exponent rj, 
z = 2 — 77 (van Hove theory) or z = 2 + crj where c = 0.7261(1 — 



1.69e ■ 



(in the e-expansion). 



length should be dropped out in the relation at the crit- 
ical point.) These two behaviors should match smoothly 
at qt; 1. The extrapolation from the hydrodynamic 
regime gives a relation such as Eq. (|2.8|) . and the ex- 
trapolation from the critical regime gives u ~ £ -z , so 
the comparison gives z. It is in this sense that we write 
t ~ £ 2 . Although we ignored this issue in the discussion 
of the model A exponent, a similar remark applies there 
as well. 

As we discuss below, an example of model B is the 
N = 4 SYM in the large- N c limit. 



III. CRITICAL BEHAVIOR OF LARGE- iV c 
GAUGE THEORIES 

A. Critical behavior of the M = 4 plasma: known 
facts 

The N = 4 SYM has the SU(4) R-symmetry, which is 
rank 3, so one can add at most 3 independent chemical 
potentials. We focus on the case of a single chemical 
potential. The theory is dual to the five-dimensional R- 
charged black hole, which is a solution of gravity coupled 
with a £7(1) gauge field and a scalar field. 

Let us summarize the known facts about the critical 
behavior of the Af = 4 SYM in the large- N c limit derived 
from the R-charged black hole. For further details, see 
Appendix lAl and original references: 

1. The R-charge density has a singular behavior 

da- 

[l9j : the R-charge susceptibility 



— 

\dfiy T 



N 2 T§ 2 + 5k-k 2 



(3.1) 



diverges at k = 2, where k is a parameter related 
to the R-charge density [see Eq. (|A4c[) ] . T is the 
temperature which is related to To, the tempera- 
ture for k = [see Eq. (|A4d|) ]. So, k — 2 is the 
critical point, and the R-charge density is the or- 
der parameter. 

2. The R-charge conductivity remains finite at the 
critical point : 



A 



(k + 2) 2 A C 2 T 



6471-^1+^ 
and the diffusion constant is given by 
1 (1 + k/2) 3 2-k 



D 



2ttT 1 + k 2 + 5k - k 2 



(3-2) 



(3-3) 



from Eqs. (|2.5|l and (|3~T|) . This implies that 
the critical slowing down indeed occurs in the 
AdS/CFT duality, i.e., D 0. It has been shown 
that z = 4 — rj with 77 = consistent with the model 
B prediction Q. 
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3. The shear viscosity remains finite at the critical 
point: 



fj = f— 2- v/1 + k 



(3.4) 



This is related to the universality of the shear vis- 
cosity fj/s = l/(47r) pol - l23j . Since the entropy 
density is the first derivative of the free energy, it 
is continuous across the second-order phase transi- 
tion. Thus, the universality implies that fj is also 
continuous across the phase transition. 

These results suggest that the M = 4 SYM belongs to 
model B rather than model H below in the large- N c limit. 
However, the system obviously has the conserved quan- 
tity T^ v . In such a case, one expects that the universality 
class becomes model H. We try to resolve this puzzle in 
this paper. But first we need some basics of model H. 



B. Model H 

The energy-momentum tensor affects the dynamic uni- 
versality class because it introduces the notion of the fluid 
motion which we ignored in Sec. Ill Dl Including this ef- 
fect modifies the dynamic universality class and the value 
of the dynamic critical exponent. 

The energy-momentum tensor contains the shear 
modes and the sound modes associated with the mo- 
mentum densities and the energy density. The scaling 
form argument in Sec. Ill Dl suggcsts that the sound modes 
lo ~ c s q oc l/£ have a weak singularity compared with 
the shear modes U) ~ —iD ri q 2 oc l/£ 2 [c s : speed of sound, 
Dfj := fj/(e + P)]. So, one needs to take only the shear 
modes into account|3 (Here, we assume that the speed 
of sound and the shear viscosity either remain constant 
or have only a weak singularity at the critical point.) 

Thus, the relevant hydrodynamic modes consist of the 
charge diffusive mode and the shear modes. The shear 
modes contain a new transport coefficient, the shear vis- 
cosity fj. These modes are known to describe model H dy- 
namics. Here, we summarize the main results for model 
H and its basic physics. For further details, see Ref. [l2[ 
and references therein. 

For model H, both A and fj become singular at the 
critical point: 



fj = £ x " ■ 



(3.5) 
(3.6) 



4 This intuitive argument ignores the couplings among hydrody- 
namic modes. In Appendix [Bj we explicitly show that the effect 
of the sound modes to the order parameter is negligible for model 
H. On the other hand, the slow dynamics of the order parameter 
can affect the sound modes. As a consequence, one has singular 
behaviors in the speed of sound and the sound attenuation (bulk 
viscosity) |24|. |25H . We will not consider this issue in this paper. 



These new dynamic exponents satisfy the following ex- 
pression: 



x\ + x n = 4 — d s 



V 



(3.7) 



For d s = 3, mode-mode coupling computations yield (as- 
suming 77 = 0) 



x x = 0.946, 



0.054 



(3.8) 



Namely, x v is a small number. The anomalous exponent 
77 is also a small number: 77 ~ for QCD and 77 = for 
the M = 4 SYM in the large- N c limit. (See also Sec. lHTCl 
for a related discussion.) Because of the constraint (13.71) . 
we eliminate x\ below in favor of x v and 77 which are 
small numbers. While we retain these small exponents, 
one may ignore them for a rough estimate. 

In order to know the values of x\ and x v separately, 
one needs mode-mode coupling computations or renor- 
malization group analysis. However, the constraint (|3.7j) 
can be understood from an intuitive argument which we 
will see below. But first let us consider the implications 
of the above result: 



3, xx 



1. 



For d_ 

limit, i.e., one cannot obtain x\ 
see the reason below. 



so there is no model B 
x v = 0. We will 



The singular behavior in conductivity occurs for 
d s < 4. This is related to the fact that the shear 
mode coupling to the charge diffusive mode is rele- 
vant in the renormalization group sense for d s < 4. 



3. From Eqs. (|3.5|) and (|3.7p . one gets 



D 



H 



For d a = 3 and 



D 



,~0, 

(Model B) 
l' 1 (Model H) 



(3.9) 



(3.10) 



Namely, the diffusion constant still vanishes at the 
critical point but with a different power of £. The 
scaling argument in Sec. Ill Dl determines the dy- 
namic critical exponent z as 



z = d s 
-iDq 2 cx £ 



2-d, 



(3.11) 

(q0 2 /Z ds+X ' 1 - 



For d s = 3 and x v ~ 0, z 



Model H does not have a model B limit. This is be- 
cause the mechanism of conduction differs for these mod- 
els. This is understood from an intuitive argument below. 
The argument also shows why the constraint (|3.7p holds. 

For model B, the conduction comes from the charge 
diffusion as we saw in Sec. Ill Dl However, for model H, 
there is an extra contribution to the conductivity from 
convection of fluids. The conductivity from the convec- 
tion is much larger than the conductivity from the charge 
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diffusion of model B (in the sense that the former diverges 
as £ — > oo). 

We introduced the external electric field E in Sec. Ill Dl 
Under E, a charged fluid will experience a mechanical 
force / appl = pE. The fluid will accelerate to velocity 
v at which viscous drag / v ; sc balances the applied force 
/vise + /appl = 0- P n this subsection, we assume p = 0, 
where ~ denotes an equilibrium value. Thus, p = 8p, 
where S denotes the deviation from the equilibrium. Sim- 
ilarly, v = Sv.) For a chunk of fluid with typical linear 
dimension L, 

f visc ~ -f W L d °- 2 , / appI ~ pEL d > . (3.12) 

Then, the induced current J is given by 

o 2 

J = pv ~ ^L 2 E . (3.13) 

n 

Thus, the conductivity Ah by convection is given by Ah ~ 
(p 2 /fj)L 2 , which diverges with L. The divergence is cut 
off at the scale L ~ £ because the fluctuation (p 2 ) = 
Tx/L ds is correlated at most this scale. Therefore, 

A H ~ *e~ ds , (3-14) 
V 

or 

^a h ~ x e- d ° ~ e- d -- v , (3.i5) 

which coincides with Eq. (|3.7[) . In Eq. p,13[) . the current 
really means J = pv = SpSv, so this effect is nonlinear 
in deviations. 

The above intuitive argument can be justified from the 
mode-mode coupling computation for r\ = which we 
reproduce in Appendix [Bj The computation also fixes 
the overall coefficient of Eq. (|3.14|) . For d s = 3, 

T 1 T 1 

Ah = tt^Xt > or D a = 7^7 ■ ( 3 - 16 ) 
birr/ £ birr/ £ 

C. Large- N c counting 

Generically, the conductivity has a model H (convec- 
tive) contribution as well as a model B (diffusive) contri- 
bution. The AdS/CFT result O is the model B con- 
tribution, and it scales as 0(N 2 ). From the boundary 
point of view, this is because there are 0(N 2 ) degrees of 
freedom which carry R-charges. 

From the bulk point of view, this comes from the 
fact that the bulk action is proportional to l/(167rG5). 
Namely, 

1 . The conductivity can be derived from the two-point 
correlator (J 1 J J ) using a Kubo formula. 

2. According to the GKP (Gubser-Klebanov- 
Polyakov)-Witten relation 0, [f|, the correlator is 
evaluated from the on-shell bulk action. 



3. The boundary current J M couples to the bulk 
Maxwell field Am- So, one needs to evaluate the 
on-shell bulk Maxwell action. 

4. The bulk Maxwell action is proportional to 
l/(167rG5) and so is the correlator!! According to 
the standard AdS/CFT dictionary, L 2 /(16nG 5 ) = 
N 2 /(8t: 2 ). 

On the other hand, the model H contribution (|3.14[) 
scales as O(l) as N c — > oo. This is because various quan- 
tities in Eq. (j3~T4")) scale as follows [Eqs. (pTTj) and (|3~4"])]: 

X,V = 0(N 2 ), £ = 0(1). (3.18) 

X and fy scale as 0{N 2 ) from similar reasons as As- The 
scaling of £ can be understood as follows. One way to 
derive the correlation length £ is to use the quasinormal 
mode method. In this case, one solves the bulk Maxwell 
equation to find a pole. At tree level, the Newton con- 
stant G5 appears only as the overall coefficient of the 
bulk action, so the bulk Maxwell equation does not con- 
tain G5, and £ cannot depend on N c . Alternatively, one 
can derive £ from the static correlator (pp) . The correla- 
tor has a structure (pp) ~ 1/(1 + q 2 £, 2 ), where the quasi- 
normal mode appears as a pole. Although the correlator 
itself is proportional to N 2 , the pole is not influenced by 
N c . 

Therefore, we conclude that the model H behavior is a 
sublcading effect at large- iV c : 

A B ~ 0(N 2 ) x e , (3.19) 
A H - 0(1) x ^-ds-v-xn . (3.20) 

The essential reason for the l/iV^-suppression is the pres- 
ence of the shear viscosity fj in Eq. (|3.16p . which is 
0(N 2 ). For model H, the effect of the convection is more 
effective than the charge diffusion since the former di- 
verges as £ — > 00. However, in large- N c gauge theories, 
the shear viscosity is larg^l and 0(N 2 ) so that the con- 
vection turns out not to be very effective. 

The anomalous exponent rj vanishes for the M = 4 
SYM in the large-iV c limit. However, static exponents 
arc also modified by l/TV^-corrections (Sec. IIVC[) . As 
seen from Eq. (|3.20p . this effect just changes the power 
of £ slightly and does not affect our main result. For 
simplicity, we ignore rj and x v in the following discussion. 



5 More generally, if the bulk Maxwell action takes the canonical 
form 

i J d 5 x^F MN , (3.17) 
x = 0(g~ 2 ), and Eq. J3T20]) changes as A H ~ 0(l/(g 2 N 2 )) X 

6 One often says that the quark-gluon plasma and gauge theories 
in the large-Af c limit has a very small viscosity, but it really 
means that fj/s or the shear viscosity for a degree of freedom is 
small. The shear viscosity itself is large due to a large number 
of degrees of freedom. 
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The R-charged black holes have been constructed for 
d s 3. When one realizes the black holes using simple 
brane systems, black holes with d s — 2,3 and 5 are par- 
ticularly important (corresponding to the M2, D3, and 
M5-branes, respectively). In the mode-mode coupling 
theory, the shear mode coupling to the charge density 
is relevant for d s < 4 if fields have the canonical scal- 
ing dimensions. Thus, the singular behavior of A may 
be invisible for the d s — 5 case even if one takes the 
l/iV 2 -effect into account. 



D. Mode-mode coupling estimate of l/JV,?-efFect 

In the AdS/CFT duality, G 5 cx 1/iV 2 . Because a su- 
pergravity 1-loop adds a weight G5, the l/A 2 -effect is 
translated into a 1-loop supergravity effect. Let us es- 
timate the model H conductivity for the AT = 4 plasma 
using the result of mode-mode coupling theory, which is 
a prediction for 1-loop supergravity. 

The mode- mode coupling result is given in Eq. (|3.1GI) . 
Reference Q estimated the correlation length numeri- 
cally. The correlator 



G(q) 



(3.21) 



1 + (qO 2 
has the pole at q 2 = — £~ 2 , where 

C 2 -0.12(2 - k){2ttT c ) 2 + O ((k -2) 2 ) . (3.22) 
Using Eqs. (|3~TB1) and (|cT2"2"]) , one gets 

Dr „ „ iiVol2^(2- K ) 1 / 2 (3.23) 

27tt 2 N 2 T 2 £ 27tt N 2 T c y ' 1 ' 

near the critical point k = 2. The corresponding conduc- 
tivity is 



32. 
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1 



A H ^0.12T 2 ^^V0T2T C(2 _ K)1/2 . 
This should be contrasted with model B results: 



Dr 



1 



6ttT c 
1 



(2 - k) 



1 



1 



24tt 3 0.12T c 3 £ 2 



A b ~ —NX 

57T 



(3.24) 



(3.25) 
(3.26) 



IV. DISCUSSION 

A. l/iV c -effects in the AdS/CFT duality 

In the AdS/CFT duality, there are various phenomena 
which arc never visible in the large- N c limit. Examples 
include "long-time power-law tails" in the relaxation phe- 
nomena 26]. the deviation from the mean-held behav- 
ior (Sec. IIVC|) . and the symmetry restoration for low- 
dimensional symmetry breaking (due to the Coleman- 
Mermin- Wagner theorem). We have seen that mode- 
mode coupling in the dynamic critical phenomena is an- 
other example of l/iV c -effects. Our result indicates that 



the dynamic universality classes of large- N c gauge the- 
ories are simpler than the usual condensed-matter sys- 
tems. 

In order to determine the dynamic universality class, 
hydrodynamic variables and conservation laws provide 
important clues. But they do not determine the dynamic 
universality class uniquely. It is also important to take 
the mode-mode coupling strength into account: this is one 
lesson of our analysis. We have seen that mode-mode 
couplings are suppressed in 1/N 2 although they are rel- 
evant operators. 

There are various attempts to explore l/iV c -effects in 
the AdS/CFT duality. See Ref. [27| for long-time power- 
law tails and Ref. |28[ for symmetry restoration. Sim- 
ilarly, it would be interesting to carry out bulk 1-loop 
computations for dynamic critical phenomena. 



B. Implications to heavy-ion collisions 

So far we focused on the M — 4 SYM in the large- N c 
limit. In this subsection, we discuss the implication to 
QCD. 

According to lattice results, the finite temperature 
transition from the hadronic phase to the quark-gluon 
plasma phase (at zero baryon chemical potential) is not 
a phase transition but rather a smooth crossover. On 
the other hand, the finite density transition at T = 
is believed to be a first-order phase transition, so there 
must be an end point of the first-order phase transition 
line somewhere in the phase diagram, which is a critical 
point. The precise location of the critical point is still 
unknown, and one of major goals of future RHIC exper- 
iments is this critical end point search. 

In real experiments, singular behaviors at the critical 
point are limited by the finite size and finite time effects 
[291 ] . The former limits £ by £ < (size) and the typical 
size of the plasma is 10 fm. However, the latter effect, 
finite evolution time is more important. Since r — £ z , £ 
is limited by £ < (time) 1 / 2 . The finite time effect is more 
important because z > 1 and the typical evolution time 
is the same order as the spatial size. 

Son and Stephanov argued that QCD belongs to model 
H at the critical end point [16| . The baryon number 
density corresponds to the R-charge densityjj For model 
H, z — 3, so the finite time effect limits £ < 2 — 3 fm 
compared with the natural value of 1 fm. 

In principle, model B and H behaviors are distinguish- 
able since diffusion constants scale with £ differently: 
Dh/Db oc £. (Here, we focus only on the baryon con- 
ductivity or baryon diffusion constant since the shear vis- 



7 Strictly speaking, QCD has an additional hydro-dynamic mode 
compared with model H, which is the QCD chiral condensate. 
However, Son and Stephanov have shown that only one linear 
combination of the baryon number density and the chiral con- 
densate becomes truly hydrodynamic. 
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cosity has only the mild singularity.) But in real exper- 
iments, one does not have enough size and time to see 
the divergence of £, and the model H effect is suppressed 
by This raises a question whether the model H 

behavior is distinguishable from the model B one in real 
experiments. 

In other words, in order to know which effect is dom- 
inant, the inspection of ^-dependence is not sufficient. 
Of course, for QCD, N% is not small but not very large 
cither. The other parameters and numerical coefficients 
would become important to compare model B and H ef- 
fects. Naively, one may imagine that numerical coeffi- 
cients are not very relevant to the comparison because the 
coupling-constant dependence would overwhelm numeri- 
cal coefficients in the strong coupling limit. One would 
expect Db,t] — > at strong coupling[f] Since Z?h oc 1/fj, 
the coupling-constant dependence would magnify the 
model H effect. However, in the AdS/CFT computa- 
tions, both Db and fj remain finite in the large- N c limit. 
This is the reason why it is necessary to reexamine the 
issue. 

It is beyond the scope of this paper to make a realistic 
estimate for QCD. Here, we compare model B and H 
effects using the AT = 4 results. From Eqs. (|3.23|) and 

(Em 



Dn 



256tt T c 
0.12-4£ 

9 Nr 



(4.1) 



Thus, the model H effect becomes dominant when 
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Nr 



256tt • 0.12 T c 



0.75 



(4.2) 



In the last expression, we used N% = 8. Thus, for the 
SU(3) TV = 4 SYM, the 1 //^-suppression does not have 
a significant effect, and the model H effect becomes com- 
patible with the model B effect as soon as £ > 1/T C . 

The above result is valid for the largc-iV c M = 4 
SYM and not for QCD, but one may be tempted to 
use the QCD transition temperature for T c , Since he ~ 
197 MeVfm and 197 MeV is not far from T c , the charac- 
teristic length scale at T c is T" 1 ~ 1 fm0 Using T^ 1 = 1 
fm, one gets £ S> 0.75 fm. It would be interesting to 
make a realistic estimate for QCD. 



For example, perturbative gauge theories give -Db > V ^ 
V(5ym m 9YM) a * zcro chemical potential [30ll , In the naive 
strong coupling limit <?ym — * 00 1 they vanish. 
Because the M = 4 SYM is scale invariant, only T c (or /j, c ) is the 
independent scale. This is not the case for QCD, and a realistic 
estimate would involve fi c as well. Also, currently the location of 
the QCD critical end point is not known very precisely. See, e.g., 
Ref. I3l| for a summary of theoretical results. For simplicity, we 
use T7 = 1 fm for the above estimate. 



C. Static critical exponents of the J\[ — 4 plasma 

The AT = 4 SYM has the following static critical expo- 
nents in the large- N c limit: 



( 2'2'2"' 2 'i'° 



(4.3) 



In this subsection, we discuss the implications of these 
exponents. 

Recall some basic facts about the static critical phe- 
nomena. In mean-field theories, the statistical fluctua- 
tions of the order parameter are ignored. The effect of 
fluctuations become more important in low spatial di- 
mensions, so mean- field exponents, which do not depend 
on the dimensionality, can be modified in low dimen- 
sions. On the other hand, the effect is less important 
in high spatial dimensions. So, mean-field exponents are 
trustable for high enough d s > d c , where d c is known as 
the upper critical dimension. 

One characteristic behavior for d s > d c is the break- 
down of the hypcrscaling relation (|2.2I) . The hyperscaling 
relation often fails above d c . This is because the so-called 
dangerously irrelevant operators may exist in free energy. 
In mean-field theories, static exponents are independent 
of d s , so the hyperscaling relation obviously fails. For 
mean-field theories, the hyperscaling relation can be sat- 
isfied only for d s = d c . 

Keeping these in mind, let us go back to the N = 4 
exponents (|4.3[) . First, these values are unconventional 
compared with the Ginzburg-Landau theory (|2.1[) . Sec- 
ond, one should regard them as an unconventional but 
mean-field result since we are taking the large- iV c limit F°l 
In the large- N c limit, fluctuations are suppressed so that 
mean-field results are exact. More explicitly, the effect 
of fluctuations can be estimated from the Ginzburg cri- 
terion. Fluctuations are suppressed when 



((p-(p)) 2 )«(p- Pc )' 



(4.4) 



Since ((p - (p)f) = 0{Nl) and (p- p c ) 2 = 0(7Y C 4 ), 
the Ginzburg criterion is always satisfied in the large- 
N c limit. As a mean-field result, one expects that the 
exponents are independent of d s . For holographic super- 
conductors, this was shown explicitly Q- The vanishing 
anomalous exponent rj = is another indication of mean- 
field results since the nonvanishing rj usually comes from 
fluctuations. 

Third, although the exponents are unconventional, 
they satisfy the static scaling relations. However, the 
hypcrscaling relation is not satisfied, which is another 
indication of a mean-field result. 



10 In this paper, we draw a distinction between mean-field theory 
and the Ginzburg-Landau theory. We use "mean-field theory" 
when the statistical fluctuations of the order parameter are ig- 
nored. The free energy of such a theory may or may not take 
the Ginzburg-Landau form. 
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Now, the story should change at finite N c . Using the 
static critical exponents, the Ginzburg criterion tells that 
fluctuations are suppressed for d > d c , where 



2/3 + 7 



(4.5) 



Using the Af — 4 exponents, d c — 6. In fact, the hy- 
perscaling relation would be satisfied at d c . Thus, the 
exponents (|4.3[) should be modified for d s = 3 at finite 
N c . The static exponents should change so that the hy- 
perscaling relation for d s = 3 is satisfied, i.e., 2 — a = 3za 
It would be interesting to compute l/iV c -effects to see the 
deviation from the "mean-field" behavior. 

The R-charged black holes have been constructed for 
d s =/= 3. Black holes with d s = 2, 3, and 5 are partic- 
ularly important. They have the same static exponents 
(a, /3, 7, S) as in the Af = 4 plasma. The exponents v and 
T] have not been computed for d s ^ 3. Suppose that the 
anomalous dimension rj vanishes even for d s ^ 3. The 
scaling relation then determines v = 1/4, so d c = 6 as 
well. Then, the static exponents should be modified for 
those systems as well. 
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Appendix A: R-charged black hole 

The Af = 4 plasma at a finite chemical potential is 
dual to the five-dimensional Einstein-Maxwcll-scalar sys- 
tcirE]: 



Sa = 



1 



167tG 5 

1 (Vm#) 2 
3 H 2 



d 5 x^/^g 



R-E H ^F* MN 



+ 2V 



where 



2V = [H 2/3 + 2H~ 1/3 



(Al) 



(A2) 



11 Capital Latin indices M, N, . . . run through bulk spacetime co- 
ordinates (t, x, u), where (t, x) = (t, x{) are the boundary coordi- 
nates and u is the AdS radial coordinate. Greek indices p,, v, . . . 
run though only the boundary coordinates. 



and L is the AdS radius. The solution is known as the 
(single) R-charge black hole (T3 - (l5| : 



■H 1 / 3 du 2 , 



dx' 



Afu 2 



ttT ^/2k(1 + h) — (dt) M 
rl 



(A3a) 
(A3b) 



A M 

n 

H = 1 + ku, (A3c) 
/ = H-(1 + k)u 2 = (1-u){1 + (1 + k)u} ,(A3d) 



where u £ [0, 1], with u — 1 corresponding to the loca- 
tion of the horizon. The parameter n is related to the 
R-charge (see below), and To is the Hawking tempera- 
ture for the neutral black hole (k = 0). Thermodynamic 
quantities are given by 



3n 2 N 2 Ti 



(1 + k) 



P 



T 



§ I 1 + K > 

n 2 N?T§ . 

2 + K 

To , 



2vT 
ttN 2 T 3 



^ v/2k(1 + ac) 



IT To 



2k 



1 



Note that T and /j, satisfy 
2ttT 




> 2 



(A4a) 
(A4b) 
(A4c) 
(A4d) 

(A4c) 
(A4f) 

(A5) 



Equation (|A5|) is invariant under k — > 4/k. This suggests 
that there are two values of k for a given pair of (T, /z) . 
One can show that the solutions with k < 2 has the lower 
free energy, so we focus on the branch. 

These thermodynamic quantities lead to the R-charge 
susceptibility and the specific heat at constant [i: 



dp 



N 2 T 2 2 



5k 



T 



df 



2 A r2rp3 (3-K)(2 + /t) 
7T iV_J n — , 

c v / TT^(2- K ) 



(A6) 
.(A7) 



They diverge at k — 2 corresponding to 7rT//i = 1, where 
the system undergoes a second-order phase transition. 

Using the conventional definition of static critical ex- 
ponents, the system has 



111 1 
- - - 2 - 

2' 2' 2' ' 4' 



(A8) 



The evaluation of v and r\ requires the correlation func- 
tion, which has been computed in Ref. 
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Appendix B: Mode-mode coupling computation of 
model H 



In this appendix, we present the mode-mode coupling 
computation for relativistic fluids. Our computation is 
similar to the one in Refs. [32l. l33j. 

In hydrodynamics, the spatial parts of conserved cur- 
rents are not independent variables; rather they are de- 
termined from conserved charges via the constitutive 
equations. For the linear deviation from the equilibrium, 
the constitutive equation for the charge current is given 
by 

Ji(t,x) = -Dd iP (t,x) . (Bl) 
At quadratic order, one can add a nonlinear term: 

Ji(t, x) = -Dd iP {t, x) + -p(t, x)Tn(t, x) , (B2) 

w 



where 



n(t,x) :=T°(t,x) , 
w := e + P . 



(B3) 
(B4) 



The quantities with ~ arc the values at a chosen equilib- 
rium state. The nonlinear term of Eq. (|B2[) represents 
the effect of convection and gives strong enhancement of 
conductivity near the critical point. 

The conductivity can be derived from the Kubo for- 
mula: 



l — f dtfdx e- ,q x (J(t, x) ■ J(0, 0))(B5) 
sT Jo J 



d 

X B + 



1 



d,w 2 T 



dt / dxe~ lqx 



x (p(t, x)ir(t, x) ■ p(0, O)tt(0, 0)) . (B6) 

The mode-mode coupling theory evaluates the four-point 
correlator as follows. First, approximate the four-point 
correlator into the product of the two-point correlators 
as 



AX{q) := X R (q) - X B 



dt / dxe~ lq - x 



d s w 2 T J 

x {p{t, x)p(0, 0)) {-Kit, x) ■ tt(0, 0)) .(B7) 
We will use the Fourier-transformed correlators 



( P (t,x)p(0,O)) 



dp 



2 -<{7r i (t,x),7r J -(0,O)}) =/ 



(27T) 

dp 



-e^ x G p {t,p) ,(B8) 
e*-G«(t,p) (B9) 



(The symmetrized correlator ({ir^x), 7Tj(0)})/2 and the 
Wightman correlator (7rj(x)7Tj(0)) make no difference for 
classical fields.) Second, use the results of the linearized 



theory for the correlators: 



G p (t,p)^e- D ^T Xp :=e' D ^ 



+ p t p j e-~ t ° p |t|/2 cc 



(Pest)} 



where p = p/p and 



In ■= — 

w 



Is 



,(B10) 
(Bll) 

(B12) 
(B13) 



The first and second terms of Eq. (|B1 1[) come from the 
transverse shear modes and the longitudinal sound mode, 
respectively. Finally, the timescale of p is much slower 
than that of 7r, so we will set D = 0. Then, (p := q — k) 



AX(q) 



1 



d*w 2 T 



dt 



dk 



T 



dt 



d s w J " J (2Tr) d - 

+e-^ p2w/2 cos(pc s t) 



——G^^Gij^p) 
\d s ~ l)e-^ 2 l t l 



■Xk 



T 



dk 



d s w J {2-rr) d - 



Xk 



4-11 



In P 2 7s P 2 + (2c s /7 s ) 2 



We henceforth perform the integration for d s = 3 and 
q = 0. One obtains 



AA 



T X 1 



6ir(w lv ) £ 67r(w 7s ) £[1 + (2c s £/ 7s )] 



(B14) 



As long as c s £,/j s ^> 1, the longitudinal sound mode 
is negligible as expected, and only the transverse shear 
modes contribute to AA: 



T 1 

Xr = X B + 7-^X7 
birr] £ 



(B15) 



A similar computation can be done for 5^0. The 
transverse shear part gives 



AA(«) 



K(qO . 



(B16) 



where 



K(x) :-- 



Ax 2 

which behaves as 



{1 + x 2 + (x 3 - x' 1 ) tan" 1 x} , (B17) 



K(x) 



(s< 1) 

_L<r- 

8 



*fx (X> 1) . 



(B18) 
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The function x 2 K(x) is the so-called Kawasaki function 
K (x). Thus, the dispersion relation u) = —iXXq 1 Q 2 be- 
comes 

uj~-i-^\q 2 (g£«l), (B19) 
birr] £ 

~ -i-^rg 3 (q£ > 1) • (B20) 



Equations (|B19[) and (|B20[) are an explicit realization of 
the scaling argument in Sec. Ill D[ In the hydrodynamic 
regime q£ <C 1, lj oc g 2 /£ ~ (<z£) 2 /£ 3 - O n the other hand, 
in the critical regime g( > 1, w a q 3 ~ (?£) 3 /£ 3 - They 
match smoothly at g£ ~ 1. 

A similar computation of the energy-momentum tensor 
correlator determines fjn and its exponent Xjj (|3.6|) . from 
which one can determine the exponent x\. 
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